Code
library(`osmdata`)
library(sf)
library(tmap)Plot wind turbines from OSM.
Load libraries
library(`osmdata`)
library(sf)
library(tmap)Get data from OpenStreetMap with osmdata, and transform to a sf object.
query <- opq(bbox = "Rogaland") |>
add_osm_feature(key = "generator:source", value = "wind")
data <- osmdata_sf(query)
wind_turbines <- data$osm_points Interactive map with tmap.
tmap_mode("view")
tm_shape(wind_turbines) +
tm_dots(col = "#0072B2")